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Abstract 

We discuss the use of a class of exact finite energy solutions to the vacuum 
source-free Maxwell equations as models for multi- and single cycle laser pulses 
in classical interaction with relativistic charged point particles. These compact 
solutions are classified in terms of their chiral content and their influence on 
particular charge configurations in space. The results of such classical interactions 
motivate a phenomenological quantum description of a propagating laser pulse in 
a medium in terms of an effective quantum Hamiltonian. 


PACS numbers: 03.50.De, 42.50.Ct, 42.55.Ah, 


1 Introduction 

Advances in laser technology have made possible the exploration of physical processes on 
unprecedented temporal and spatial scales. They have also opened up new possibilities 
for accelerating charged particles using laser-matter interactions. Multi- and single cycle 
high intensity (10 10 — 10 15 Wcm -2 ) laser pulses can be produced using Q-switching or 
mode-locking techniques [I]. Pulses of even higher intensity ( 10 21 Wcm 2 ) could 

accelerate charged particles such as electrons to relativistic speeds where radiation 
reaction and quantum effects may influence their dynamics |2|. Lower intensity pulses 
have also been used as diagnostic tools for exploring the structure of plasmas in various 
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states |3,4|. In order to interpret experimental data involving classical laser interactions 
with both charged and neutral matter, theoretical models |5-8| rely crucially on 
parameterizations of the electromagnetic fields in laser pulses, particularly in situations 
where traditional formulations using monochromatic or paraxial-beam approximations 
have limitations (9 11 . Such theoretical models may not be compact in all spatial 


dimensions and the role of a laser pulse as a classical probe is further limited by the 
scales that it is designed to resolve. Pulse shape design characteristics are often guided 
by simulations of laser-matter interactions which incorporate the known laws of physics 
of relevance at such scales. For example, the intense experimental activity currently 
exploring the electromagnetic properties of single-cycle laser pulses with nanoscale 
objects (such as dielectric and plasmonic nanoparticles) demands efficient modelling 
tools that accommodate the spatial compactness of such pulses. Such tools will 
eventually require incorporation of quantum effects associated with these interactions 
in order to properly describe observations and yield practical applications involving 
nanostructured dielectrics and plasmonic metals 12 . Furthermore an effective quantum 


description of the laser pulse itself would offer a new simulation tool for designing more 
accurate methods of encoding quantum information. 

The quanta associated with plane time-harmonic electromagnetic fields in vacuo provide 
ideal discrete two-state systems (photons) that are routinely used as controllable qubits 
in information science. Three-level bosonic quantum systems composed of two photons 
in the same spatial and temporal configuration have also been contemplated jl31 in 
attempts to construct more efficient quantum gates for quantum communication. If 
one regards a free propagating classical single-cycle (therefore non-time-harmonic) laser 
pulse in vacuo as a spatially compact classical electromagnetic configuration with definite 
energy £ and temporal width to, one expects that when £t 0 < h its dynamical evolution, 
in both vacuo and material media, should be controlled by an effective quantum 
Hamiltonian, rather than a classical one as is done in 114 . Such quantized collective 


states could then be entangled with other quantized pulses or free photon states and 
their interaction with classical or quantized states of electrically neutral continua (e.g. 
optically inhomogeneous and anisotropic dielectrics or plasmas) or charged matter (e.g. 
trapped ions [l3|) may be worthy of investigation for technological applications such as 
quantum computing and encryption. 

In this article we first discuss a viable methodology for parameterizing a particular class 
of propagating solutions to the source-free classical Maxwell equations in vacuo that 
offers an efficient means to explore the classical effects of compact laser pulses on free 
electrons in dynamical regimes where quantum effects are absent. The parameterization 
is constructed from a remarkable class of explicit solutions of the scalar wave equation 
found by Ziolkowski 15 -19 following pioneering work by Brittingham 20 and Synge 


21 . Solutions in this class are parameterized in terms of three real constants that 


are sufficient to completely determine the characteristics of any freely propagating 
laser pulse in full accord with Maxwell’s equations in free space. They describe 
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solutions with finite total electromagnetic energy, electromagnetic fields bounded in 
all three spatial directions and experimentally distinguishable chiral configurations. 
With simple analytic structures their diffractive properties can be readily calculated 
together with the behaviour of relativistic charged particle-pulse interactions over a 
broad parameter range without recourse to expensive numerical computation. Based 
on this behaviour and by analogy with the effective Hamiltonian theory of diatomic 
molecules, we are led to construct a Hilbert space on which to describe certain quantum 
states of an electromagnetic pulse undergoing unitary evolution generated by an effective 
phenomenological quantum Hamiltonian in a medium. Such a Hamiltonian, defined in 
terms of a set of parameters associated with its medium interaction, serves as an effective 
quantum model for the quantum pulse evolution. 


2 Parameterizing Compact Electromagnetic Laser 
Pulses 

If a complex scalar field a satisfies □ a = 0 and is any covariantly constant (degree 
2) anti-symmetric tensor field on spacetime (i.e. = 0 for all /q u, 8 — 0,1, 2, 3), then 

the complex tensor field = d^A v — d u A fl satisfies the source-free Maxwell equations 
in vacuo with: 


A v = ^(an^e^VH (!) 

where \g\ is the modulus of the determinant of the spacetime metric tensor field g with 
components g )iu and e 1 ^ v denotes the Levi-Civita alternating symbol. In the following, 
we restrict to Minkowski spacetime, in which case the components YL^ U can be used to 
encode three independent Hertz vector fields and their dual^j] 

General solutions to □ a = 0 can be constructed by Fourier analysis. In cylindrical polar 
Minkowski coordinates {t,r, z,6}, axially symmetric solutions propagating along the z- 
axis have, for z > 0, the double integral representation a(t, r, z) = du e~ luJt a(co, r, z ) 
where: 

a(w,r,z) = J'kUk)Mkr)e^W^ dk + J kf u (k) J 0 (fcr)e"V fc2 -(“) 2 dk 

in terms of the zero order Bessel function and the speed of light in vacuo c. 

Conditions on the Fourier amplitudes f w (k) can be given so that the Hertz procedure 
above gives rise to real singularity free electromagnetic fields with finite total 
electromagnetic energy. A particularly simple class of pulses that can be generated 

Hn the language of differential forms on Minkowski spacetime A = *d(all), F = dA where 
d * d a = 0, the 2-form II satisfies VII = 0 and * denotes the Hodge map associated with g. 
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Figure 1: Power profile P of the (CM, 1) laser pulse at Z = 0, T = 0 with parameters 
{A = 600, tfi = 1, 2 = 1000, <f> = 0.001, E = 1} 


in this way follows from the complex axi-symmetric scalar solution: 


a(t,r,z ) = 


r 2 + (0i + i(z - ct )) (0 2 - i(z + ct )) 


( 2 ) 


where 0 ), 0 i ,02 are strictly positive (real) parameters with physical dimensions of 
length. The relative sizes of 0i and "02 determine both the direction of propagation 
along the z—axis of the dominant maximum of the pulse profile. When 0i 02, 
the dominant maximum propagates along the z—axis to the right. The parameter £o 
determines the magnitude of such a maximum. The structure of such solutions has been 
extensively studied in 22,23 in conjunction with particular choices of fl „ v together with 


generalizations discussed in 24,25 


In general the six anti-symmetric tensors with components dCW in a Minkowski 
Cartesian coordinate system are covariantly constant and can be used to construct 
a complex eigen-basis of antisymmetric chiral tensors II s ’ K , with s G {CE,CM} and 
k G { — 1, 0,1}, satisfying 


O z hP’ K = kIP’* * (3) 

where the operator O z represents 6 rotations about the axis generated by —idg 
on tensor^] These in turn can be used to construct a complex basis of chiral eigen- 
Maxwell tensor fields . The index s indicates that the CE (CM) chiral family 
contain electric (magnetic) fields that are orthogonal to the z—axis when k = 0. 
The chiral eigen-fields F^° inherit the axial symmetry of a(t,r,z ) while those with 

Tn terms of the Lie derivative, O z = —iCd e and n CE ±1 = d{x ± iy) A dt , n CE ’° = dz A dt , n CM,K = 

*n CE > K where x = rcos(d),y = rsin(0) 
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Figure 2: Three-dimensional spacecurves for particles subject to an incident (CM, 1) 
laser pulse (left), (CM,0) laser pulse (centre) and (CM, —1) laser pulse (right) with 
parameters {A = 600,dR = 1,T 2 = 1000,$ = 0.001, H = 1}. Each particle has initial 
velocity (-R(O) = 0, 0(0) = 0, Z( 0) = -A,}. The shaded circular disc region indicates the 
initial spot size (R = 10000 for (CM, ±1) laser pulses and R = 20000 for a (CM, 0) laser 
pulse) relative to the black markers on the spacecurves that denote the initial positions 
of the charged test particles. 

k = ±1 do not. The directions of electric and magnetic fields in any of these Maxwell 
solutions depend on their location in the pulse and the concept of a pulse polarisation 
is not strictly applicable. The chiral content as defined here can be used in its 
place. Non-chiral (complex) pulse conhgurations can be constructed by superposition: 
F lJV = Y, s C S,K with arbitrary constant complex coefficients C S,K . 

The energy, linear and angular momentum of the pulse in vacuo can be calculated from 
the components T^ u of the Maxwell stress-energy tensor T^ u = — \fj IJ ,uF o ^F a 0 ~ F^ a F a u 
where T jiv = R e(F^ u ). If e and b denote time-dependent real electric and magnetic 
3-vector fields associated with any pulse solution F^, its total electromagnetic energy 
J, for a fixed set of parameters and any z, is calculated from 

J = — [ dt [ (e x b) • dS (4) 

Po i-oo Js 

where S can be any plane with constant z = Zo > 0. For spatially compact pulse fields 
in vacuo this coincides with the total pulse electromagnetic energy 

p poo p2n poo 

£ = pdV= dz dO j rdr p(t , r, z, 9) (5) 

JV J~oo Jo Jo 

where p= \ ^eo e ■ e + ^ j is integrated over all space V. This follows since V-(exb) = 
— po d t p. To correlate J with other laser pulse properties and the choice of parameters, 
we bring the pulses F for various values of s and k into classical interaction with 
one or more charged point particles. The world-line in spacetime of a single particle, 
parameterized in arbitrary coordinates as = £ m (t) with a parameter r, is taken as a 
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solution of the coupled non-linear differential equations 


Mr) = —^Mr))V(T) ( 6 ) 

m 0 c z 

in terms of the particle charge q and rest mass mo, for some initial conditions £(0), V (0), 
where the particle 4-velocity satisfies V v V v = —1 and its 4-acceleration is expressed in 
terms of the Christoffel symbols as A^r) = d T V^(r) + Vs(r) Vp(r) r < 5 / 5 /i (^(r)). In 
the following, radiation reaction and inter-particle forces are assumed negligible. From 
the solution £(t) one can determine the increase (or decrease) in the relativistic kinetic 
energy transferred from the electromagnetic pulse to any particle and the nature of its 
trajectory in the laboratory frame for different choices of s and k. 


3 Interactions of Compact Laser Pulses with Matter 

The analysis of the previous section can be used to investigate the behaviour of a classical 
laser pulse interacting with electrically charged matter and motivates a model for a 
quantum laser pulse interacting with electrically neutral matter. 

In order to facilitate the classical behaviour, we reduce the above equations of motion to 
dimensionless form and fix the physical dimensions of the fields involved. The Minkowski 
metric tensor field g = g^dx 11 dx u (with g liu = diag(—1,1,1,1)) in inertial coordinates 
x° = ct,x 1 = x,x 2 = y,x 3 = z) has SI physical dimensions [ L ] 2 . The SI dimension 
of electromagnetic quantities follows by assigning to eo F^dx^dx 1 ^ in any coordinate 
system the physical dimension of charge. Furthermore, in terms of Minkowski polar 
coordinates {t,r, z,9}, introduce (for ease of graphical visualization) the dimensionless 
coordinates {/? = r/<M 0 , T = ct/£ 0 , Z = ^/S£ 0 } and real dimensionless parameters 
A, Tj = ipj/i o (j = 1,2) where [dt,-] = [<f>] = [S] = 1, [£ 0 ] = [L\. Then with the 
dimensionless complex scalar held a(T,R,Z) = a(t,r,z ) and Greek indices ranging 
over {T, R, Z, 6} with e T ’ R ’ z ’ e = 1 , we write 

4s = AAa ? (sn /w ) G'45 (7) 

for a choice of dimensionless covariantly constant tensor II^(T, R, Z, 9) = II ^(t, r, z, 9), 
so that [eo dx IJ ] has the physical dimension of electric charge. The total power density 
P and total energy density S are now defined by 

/ OO roo 

dT dR d9P(T,R,Z,9) 

-oo JO Jo 

/ OO 

dZS(T,Z). 
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The real parameter A controls the strength of the electric and magnetic fields in Ffifi 
for fixed values of the real parameters Ti, T2, <E>, S and the overall scale £0 will be 
fixed in terms of the total electromagnetic energy of the pulse. For a choice of such 
parameters the associated real fields e and b enable one to calculate a numerical 
value T such that J = ZqT. The diffraction of the pulse peak along the z—axis 
can be used to define a pulse range relative to the maximum of the pulse peak at 
z = 0. To this end, the density £(T,Z) defines the dimensionless range Z rg by 
£(0,0 )/S(Ti, Z rg ) = 2, where the peak at Z — Z rg > 0 and T = Ti > 0 is half the 
height of the peak at Z = 0,T = 0. If during the interval [0, T\] the pulse propagates 
with negligible deformation in Z, one may estimate its width Z w at half height and 
the dimensionless pulse axial speed f3 = Z rg /T\. This yields the dimensionless pulse 
duration or temporal width T 0 = Z w /f3. From these dimensionless values one deduces 
the pulse SI characteristics in terms of £ 0 and hence J. If the picosecond is used 
as a unit of time, the pulse duration becomes t 0 = £qTq/c = £ 0 Z w /(f3c) = iV10~ 12 
s for some value N and hence £ 0 = (c/3N/Z w ) 10~ 12 m, J = (TficN / Z w ) 10~ 12 J, 
z w = ’Ecfin 10 -12 m and z rg = £oEZ rg = (E(3cNZ rg /Z w ) 10~ 12 m. A dimensionless 
spot-size of the pulse at Z = Z 0 > 0, T = Z 0 /(3 is then determined by the behaviour 
of P(R, Zq/ f3, Z 0 , 6). At each value of Z 0 this function of R and 0 has a clearly defined 
principal maximum. If one associates a circle of dimensionless radius R S (Z 0 ) with 
such a maximum locus it can be used to define a spot-size at z = z 0 with radius 
r s (zo) =£o&Rs{Z 0 ) = (c/3N$R s (Z 0 )/Z w )10~ 12 nr. Figure [l] displays a clearly pronounced 
principle maximum in the power density profile P as a function of X = Rcos(6) and 
Y = R sin(0) at Z = 0, T = 0 for a specific choice of the parameters (A, Ti, T2, $, 5). 
The same parameter set is used to numerically solve ([6]) for a collection of trajectories 
for charged particles, each arranged initially around the circumference of a circle in a 
plane orthogonal to the propagation axis of incident CM type laser pulses with different 
chirality k. The resulting space curves in 3-dimensions, displayed in figure [2| clearly 
exhibit the different characteristic responses to CM pulses with distinct chirality values. 
The instantaneous specific relativistic kinetic energy of a particle with laboratory speed 


s 

K 

£(J) 

A 

Zrg (m) 

t 0 (ps) 

Z w (ill) 

J's(O) (m) 

I (W cm 2 ) 

CE 

0 

8.07 

40000 

1.44 

3 

8.99 x 10 -4 

0.018 

2.646 x 10 12 

CE 

±1 

7.69 

1 

2.43 

3 

8.99 x 10" 4 

0.009 

1.009 x 10 12 

CM 

0 

8.41 

1 

20 

0.899 

2 

6.00 x 10 -4 

0.012 

9.304 x 10 11 

CM 

±1 

7.69 

1 

1.80 

3 

8.99 x 10“ 4 

0.009 

1.009 x 10 12 

CM 

±1 

76903 

100 

1.80 

3 

8.99 x 10 -4 

0.009 

1.009 x 10 16 


Table 1: Table showing SI laser characteristics for various (s, k) laser pulse configurations 
with parameters {db = 1, T 2 = 1000, $ = 0.001, H = 1}. 
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Figure 3: Specific kinetic energy transfer 7 — 1 to any charged test particle. On the 
left, the (CM,1) pulse has parameters {$ 1 = 1, ^/ 2 = 1000,$ = 0.001, E = 1}. On the 
right, relative values of 7 — 1 are displayed for various (s, k) pulses with parameters 
{A = l,$i = 1, $2 = 1000,$ = 0.001, H = 1}. The differences between the 
energy transfers for some pulses appear indistinguishable relative to others owing to 
the logarithmic scales employed. I 11 all cases, the charged particle has initial position 
(-R(O) = 1, 0(0) = |, Z( 0) = 1} and initial velocity (-R(O) = 0, Z(Qt) = Wp 0(0) = 0}. 


this 


v is 7 — 1 in terms of the Lorentz factor 7 given by 7 1 = \J 1 — In figure 
quantity is displayed as a function of T on the left for a charged particle accelerated by a 
fixed chirality (CM, —1) type pulse where the pulse energy is varied by changing A. On 
the right the energy transfer dependence on pulse chirality for both CE and CM type 
pulses with fixed laser energy is displayed. We deduce that the pulse momentum and 
angular momentum 26 in the propagation direction can transfer an impulsive force 
and torque respectively to charges lying in an orthogonal plane. More generally, the 
classical configurations of a high energy pulse labelled CE or CM could be distinguished 
experimentally by their interaction with different arrangements of charged matter. 
Furthermore, by a suitable choice of parameters, (CE, k) type modes can be constructed 
that yield the same physical properties (Jz rg , z w , /3) for all k. Similarly the (CM,k) 
type modes yield a n independent set with physical properties distinct from those 
determined by the (CE,k) modes. The pulse group speed magnitudes (as defined 
above) of all these configurations are determined numerically and are bounded above 
by the value c. To illustrate some of these statements, Table [3] summarizes the SI laser 
pulse characteristics determined by solving the system of equations in (| 6 ]) for a specific 
choice of {$ 1 , $ 2 , $, 5} and various values of A, selected to demonstrate the physical 
characteristics of existing laser pulses. We conclude that both the CE and CM type 
pulses can propagate localized packets with linear and angular momentum analogous 
to the quantum characteristics of propagating wave packets describing particles with 
different masses. 

These observations suggest that, for electromagnetic micropulses with £t 0 < h, 
a parameterized effective quantum Hamiltonian H describing a particular pair of 
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interacting massive point particles may provide the simplest effective description of 
general non-stationary quantum states of a laser pulse. Such states are then defined 
as elements of a complex Hilbert space T~L = X 2 (R 3 , C 3 ) 0 X 2 (R 3 , C 3 ) where each factor 
denotes the space of complex square integrable 3-component vectors on R 3 carrying 
irreducible finite dimensional representations of the 3-dimensional rotation groupj^] 
Such a construction is analogous to the Hilbert space "H vi b 0 'Hmt used in Bohm and 
Loewe’s dynamic model of a vibrating and rotating diatomic molecule after removing 
its translational degrees of freedom 27 . The space "H v ib is spanned by eigenstates of 


the Hamiltonian H v ^ for a one-dimensional simple harmonic oscillator parameterized 
by a fundamental frequency u> 0 and the space H ro t is spanned by eigenstates of the 
Hamiltonian H vot for a one-dimensional rigid rotator parameterized by a constant 
moment of inertia parameter X. The simplest model for the diatomic molecule takes for 
its Hamiltonian: 


H mo \ — H v ib 0 / + / 0 H mt + tylTvib 0 Wbot 

in terms of the three phenomenological real parameters Q and where / denotes the 
identity operator on the relevant space. For certain diatomic molecules this idealized 
model (and its refinements) accounts for well-established empirical formulae for their 
electronic energy states. States in the space % must describe the translational modes of 
a quantum packet as well as possible mode-mixing. To this end we model their unitary 
evolution generated by an effective Hermitian Hamiltonian of the form: 

H = XX 0 (xi , x 2 ) + 77 int (xi, x 2 ) 


where 

with real parameters i^i, > 0. For quantum pulses that are deemed relativistic, one 
can include relativistic corrections by the replacement 

—► yf-tPc 2 ^ + /i 2 c 4 - fi n c 2 

for n — 1,2 and working in a momentum representation. The structure of i7i nt (xi, x 2 ) 
depends upon the physical properties of the medium through which the pulse propagates. 
A simple model that accommodates quantum birefringence , medium anisotropy and 
inhomogeneity may be constructed in terms of parameterized Hermitian operators 
H^xO, W 2 (x 2 ): 


Hint = VFl(xi)0fH 2 (x 2 ). 

§The inner product on T~L = Hi <g> H .2 is defined by (cfi 1g> ® <t>'%) = (fi{4>2A'flUi for 

<j>n, 4>' n £ T~L n , n = 1, 2 in terms of the inner products <t>' n )'Hn on Tin 
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The six classical chiral states labelled (s, k) that evolve according to the classical Maxwell 
equations are now replaced by elements \Eq G H satisfying the Schrodinger equation: 


h d T 

= 

i dt 




with \Eq o prescribed at any time t 0 and satisfying (\Eq o ,\l/ to ) = 1. In a direct product 
matrix representation, decomposable states in % can be written ct(t, x x ) <g) /3(t, x 2 ) with 


a(t,xi) = 


^ Q!l(t,Xi) ^ 
a 2 (t,xi) 

ya 3 (i,x i) 


and /3(t, x 2 ) = 


^ Pi (t,x 2 ) ^ 

P 2 (t,X 2 ) 

y^ 3 (t,x 2 ) y 


By analogy with the description of non-relativistic spin—| qubit states, call ol the 
left qutrit component of such a bi-qutrit state and (3 the corresponding right qutrit 
component. The dynamic description of a general state ipt £ requires an explicit 
interaction Hamiltonian operator. To model interactions with a classical medium we 
envisage here a smooth fabricated meta-material with specified inhomogeneous and 
anisotropic characteristics. This is similar to the specification of the magnetic moment 
interaction of the qubit states (due to an unpaired bound s-wave electron) of a silver atom 
with a classical (inhomogeneous) static magnetic field in the Stern-Gerlach experimental 
arrangement. An important practical distinction however arises since the laser states, 
unlike electrons, are electrically and magnetically neutral and hence their interaction 
with electromagnetically neutral properties of a fabricated meta-material may be more 
versatile than in the Stern-Gerlach situation 


Consider, for example, the evolution of the left qutrit component a(t,xi) under the 
Hamiltonian: 


K 2 2 
-V 2 

2/ii X1 


I + GiS- JV(t, Xl )®7 


where the constant (3 x 3) Hermitian matrix-valued vector S satisfies the commutation 
relation S x S = iS, the 3-vector N(t, x x ) is a classical field of unit-vectors on some 
fabricated meta-material and G± is a parameter with the physical dimensions of energy. 
Such a Hamiltonian assigns a single preferred field of spatial directions N(t, x x ) for the 
interaction of quantum states with the medium. More complex media could involve 
multiple anisotropies described by a multi-directional set of unit-vector fields. 

Given a(0,X!) with 



aq(0,xi)| 2 d 3 xi 


1 
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(assuming an unbounded medium id), one solves the system: 


hdcx{t,^ci) 
i dt 


2 /ii 


V^a^x i) + Gi S ■ N(t, xi) a(t, Xi) 


for the left qutrit component at £ > 0. Multi-component wave-packet solutions can be 
constructed by the standard Fourier transform methods used to construct wave-packets 
for scalar fields. Thus in cylindrical polar coordinates (r, </>, z) defined with respect to a 
preferred time-independent, fixed z—direction N = (0,0,1) in a uniformly homogeneous 
medium, wave-packet solutions to the above equation take the form: 


x /»C 

aj(t,r,(j),z) = / 

^ J — C 


dk A 


3,™ 


( k,s ) / \/sr ds J m (rs) exp{ikz + im(j)\ Tj(t,k, s) 

Jo 


for j = 1, 2, 3 where 


Ti(t, k,s) 
T 2<yt,k,s) 
F 3 {t,k,s) 


exp 

exp 

exp 


f (|> 2 + s2 > + G ‘ 


and the amplitudes Aj }m (k,s ) are determined from the normalised initial conditions 
aq(0,r, (f, z) by Fourier and Fourier-Bessel inversion. In these expressions contributions 
to ay from terms in the sums that depend on the integer m indicate those from 


eigenstates of the orbital angular momentum operator -Cg 


29 


By contrast, terms 


in a with different j refer to spin (qutrit) contributions to ift- With the above 
Hamiltonian, stationary qutrit energy eigenstates would be non-degenerate with energy 
shifts (0, ±G*i). Generating superpositions involving left and right qutrit states in "H 
that cannot be reduced to decomposable states by a change of basis may offer a means 
to isolate and thereby control non-stationary entangled qutrits using appropriately 
fabricated meta-materials. 

In addition to such laser state “measurement” interactions with a classical medium, one 
may also include interactions with atomic quantum states. These and other uses of such 
effective quantum Hamiltonians involving the dynamics of laser states will be discussed 
elsewhere. 
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4 Conclusions 


We have constructed a basis of classical chiral solutions of the source-free vacuum 
Maxwell held equations from a simple particular solution to the complex scalar wave 
equation in spacetime and a set of covariantly constant antisymmetric tensor fields. 
Such solutions offer a simple 3-parameter description of a finite-energy laser pulse 
that provide a more accurate simulation tool for analysing laser-matter interactions 
in realistic 3-dimensional situations where plane-fronted paraxial approximations are 
inadequate. The analytic structure of such solutions enables one to readily extract 
all the standard diffractive characteristics associated with a laser pulse in free space. 
Using the classical relativistic Lorentz-force equation of motion we have also analysed 
numerically the interaction of such compact pulse solutions with charged point particles. 
This has explicitly demonstrated how laser configurations with definite chirality and 
mode-type transfer angular momentum and energy to the charges as a result of the 
interaction. From these numerical investigations we have proposed a particular effective 
quantum model for systems where the classical pulse energy S and the pulse duration t 0 
satisfy £to < h. By analogy with the effective modelling of rotating-vibrating diatomic 
molecules we have proposed a simple phenomenological Hamiltonian that may be used 
to describe quantum laser packets in free space and material media, ft is suggested that 
this Hamiltonian may have utility for simulating a novel transfer of quantum information 
and for constructing models of rapid single-cycle laser pulses interacting with quantum 
matter and classical fabricated materials containing structures below the nano-scale. 
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